Asymmetric bubble disconnection: 
persistent vibration evolves into smooth contact 
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The disconnection of an underwater bubble iUustrates how slight initial asymmetries can prevent 
the formation of a finite-time singularity. Creating a singularity by focusing a finite amount of 
energy dynamically into a vanishingly small amount of material requires that the initial condition 
be perfectly symmetric. In reality, imperfections are always present. We show a slight azimuthal 
asymmetry in the initial shape of the bubble neck excites vibrations that persist over time. As a 
result, the focusing singularity is generically pre-empted by a smooth contact. 
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Whenever a diver plunges into a swimming pool or 
a waterfall hits a river, air is entrained, often as large 
cavities that subsequently break-up into many pieces [Ij. 
Here we focus on a simple version of this disintegration 
process: the disconnection of an air bubble from a nozzle 
while it is submerged under water (Fig. la). Previously 
it was believed that the bubble neck disconnects at a sin- 
gle point via an axisymmetric implosion that ends in a 
focusing singularity glElllISlEllZlHll]. In this sce- 
nario, water outside the bubble rushes inwards to con- 
strict the neck of air, a motion converting the poten- 
tial energy during the initial instant into kinetic energy. 
Since the amount of water rushing inwards decreases to 
as the neck radius goes to 0, all the kinetic energy is 
concentrated into a vanishingly small mass at the mo- 
ment of break-up [10]. Analogous focusing singularities 
arise in models of sonoluminiscence [11 , supernova and 
shock- wave implosion [12]. 

Nearly all the previous studies analyzed the implosion 
by assuming that the focusing singularity controls the 
final dynamics. However, Keim et al. found in experi- 
ments that a slight, azimuthal asymmetry in the shape 
of the bubble neck completely transforms the dynam- 
ics [6^. An example of the different outcome is given 
in Fig. lb. Recently, Schmidt et al. [10 showed that 
the initial asymmetry transforms the final disconnection 
by exciting vibrations in the cross-section shape of the 
bubble neck. Once excited, the vibrations persist with 
the same amplitudes. Since the average radius of the 
bubble neck is decreasing to 0, the interface always be- 
comes strongly distorted as break-up approaches. These 
constant-amplitude vibrations encode a detailed memory 
of the initial state. Such a detailed memory is surpris- 
ing and important because the severely nonlinear evo- 
lution towards a singularity is commonly thought to be 
governed by convergence onto a universal dynamics, one 
independent of boundary and initial conditions and thus 
having little memory of the initial state [131 [141 [15] . More 
generally, there are evidences supporting the view that 
such memory-encoding vibrations is a common feature of 





FIG. 1: Bubble disconnection dynamics, (a) Experiment: an 
air bubble (dark area) is submerged under water and released 
from a nozzle. Bright spots are optical artifacts, (b) When 
the bubble is released by a sudden burst of gas pressure from a 
slot, the initial ribbon- like neck disconnects via a "puncture", 
instead of breaking at a point. There is a circular "hole" of 
water (light area) inside the ribbon-shaped bubble neck (dark 
area). A satellite bubble (dark dot) is visible in the water- 
filled hole. See [6] for details, (c) Model: A cross-section of 
the neck (bubble region) contracts inwards due to a radially- 
symmetric influx of water from the far- field in the exterior. 



focusing singularities [IOl[Tl[T7l[l8llia|20l|2ll|22] 

Here we focus on consequences of memory which can- 
not be addressed within the assumption of weak distor- 
tion employed by Schmidt et al [10] . We simulate the in- 
terface evolution in the final moment of break-up, when 
the surface is strongly distorted. We find that the bub- 
ble neck generically evolves towards a finite-time contact. 
Initially distant portions of the interface osculate, touch- 
ing each other at a point and having a common tangent 
at the point of contact (Fig. 2b). We also show that 
the finite-time contact can be viewed as a natural conse- 
quence of the vibrational amplitudes persisting with the 
same size while the average radius decreases to 0. 

Our analysis starts from the simplest situation: the 
break-up of a long and slender bubble neck. Previous 
works have shown that, in this asymptotic limit, every 
cross-section of the bubble neck contracts solely due to 
the influx of water in the same horizontal plane, with no 
influence from the dynamics at other heights [2 , 3] l9j . We 
therefore focus on the cross-section at the height where 
the break-up will first occur. The evolution of the cross- 
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section at a fixed vertical height corresponds to a two- 
dimensional (2D) version of the classic Rayleigh-Plesset 
collapse |23] (Fig. Ic). Stacking the 2D evolutions at 
the different heights vertically then generates the time- 
evolution of the entire neck. We also neglect the effects 
of air flow within the bubble neck and surface tension 
at the interface. Both effects can modify but do not 
change the essential features of our results. Schmidt et al. 
have shown that the 2D model described above quantita- 
tively reproduces the measured cylindrically-symmetric 
dynamics until the neck is a few microns across [10 . 

We use a polar coordinate system whose origin coin- 
cides with the center of the neck cross-section. The air- 
water interface is given by r = S{6^t). Since viscous 
effects are negligible, the exterior velocity field is irrota- 
tional, i.e. u(x, t) = V<^> where <l> is a velocity potential. 
Since the exterior flow is also incompressible V-u = 0, the 
velocity potential $ satisfies Laplace's equation V^<l> = 0. 
At t = 0, a distorted bubble- neck, modelled as a region 
of constant and uniform pressure, is immersed in water. 
To simplify the analysis, we choose the reference pres- 
sure level such that the bubble pressure p is 0. To drive 
the break-up of the bubble, we require that the exterior 
flow far from the bubble approach the form (— Q/r)e^ as 
r ^ oo. To simplify the calculations, we require that, 
27r(5, the volume flux of the far- field in- flow, remains 
constant over time. This is equivalent to requiring that 
R{t), the average radius of the bubble cross-section, de- 
creases as i^o\/(^* —t)/t^^ where Rq is the initial value 
and = Rq/2Q is the time when the average size of the 
bubble goes to 0. This square-root decrease reproduces 
the leading-order behavior of the measured break-up dy- 
namics O m [71 El ini Ho]. Schmidt has shown that the 
higher-order, slow variation in Q does not affect behav- 
ior of the vibrations at leading-order [22 . Two bound- 
ary conditions dictate the time-evolution of the velocity 
potential ^{r^O^t) and the interface S{0^t). First, at the 
rapidly accelerating surface, the normal stress exerted by 
the exterior flow must equal the bubble pressure. This 
unsteady version of the Bernoulli condition has the form 



[d^/dt+{l/2)\V'^f]\s = 0. 



(1) 



Second, the position x of a material point on the surface 
is advected by the exterior flow, 



dx 

dt 
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dt 



V$-V x = V$|s 



(2) 



We use an approach developed by Dyachenko et al. to 
solve for the interface evolution [24l [2 5 . First, the solu- 
tion of V^<l> = in the exterior is simplified by mapping 
the exterior of S{0, t) conformally onto the exterior of a 
unit circle in the complex plane w (see e.g. [26 ). A point 
on the iL'-plane is related to the location (r, 0) on the 2D 
plane by z{w,t) = re*^. The velocity potential in the 
physical plane is then given by the real part of a complex 



velocity potential ^. Second, instead of solving directly 
for the interface z = S{0^ t)e^^ and the complex potential 
^, we work with the new variables 

n{w, t) = l/{wd^z) V{w, t) = (S^^)/(5^z). (3) 

The variable V corresponds to the speed of the fluid mo- 
tion on the interface. The variable IZ does not have a 
straightforward physical interpretation, though it clearly 
is a measure of how distorted the void shape in the real 
space has become relative to the unit circle on the w- 
plane. Equations (1) and (2) now assume the form 

dtU = ^(5^7^)^{Re[7^V*]} - ^7^5^^{Re[7^V*]} (4) 
dtV = ^(S^V)^{Re[7^V*]} - ^7^S^^{|Vp} /2 (5) 

where the * symbol denotes complex conjugation. The 
integral operator A is the Cauchy integral (see e.g. 
[26j J7 ): In the simulation, we represent the initial state, 
i.e. the surface shape S and the velocity potential ^, in 
terms of IZ and V, as an expansion with N total modes, 
^ = En=o Cn/^"+' and V = En=o ^n/^"+^ The coef- 
ficients Cn and dn are computed from ^ and S. We then 
update the values of 1Z and V using (4[5). The new neck 
shape S and velocity potential ^ are in turn computed by 
numerically integrating ([3|. This formulation speeds up 
the numerical computation because the right hand side 
of equations (4j5) can be computed in TV log steps us- 
ing the Fast Fourier Transform. Results presented below 
have been obtained with A^ = 256. Using A^ = 1024 or 
more points yields no significant changes. 

Fig. 2a and 2b give the outcome of a typical simulation. 
The average radius R decreases as ^/(t^—t)/t^ (Fig. 2a). 
As the bubble region shrinks, an initially slight distortion 
becomes more and more noticeable. To illustrate how the 
distortion develops we have chosen 6 moments (labelled 
in Fig. 2a) and rescaled them by R{t). At t = the void 
shape is a nearly circular oval, given by S{0^t = 0) = 
Rq + ^2 cos 1^2 cos (2^) where Rq is the initial value of 
R, A2 the size of the azimuthal distortion and Q2 the 
initial phase. During the early moments, A2 <C R{t), 
so the cross-section vibrates, in quantitative agreement 
with linear stability results [10 , 

S{0,t) = R{t) ^ A2Cos[(l)2{R)]cos{20) (6) 

MR) = HRo/R{t))^n2- (7) 

The amplitude of the vibration A2 remains constant 
while the phase of the vibration ^2 oscillates. In Fig. 2b, 
the cross-section is initially elongated along the x-axis, 
then along the ?/-axis, and later is elongated along the 
X-axis again. As the break-up proceeds, the interface 
becomes strongly distorted {A2 ~ R{t))^ evolving out 
of the the linear stability regime. The cross-section de- 
forms into a dumb-bell shape. As time goes on, the waist 
of the dumb-bell narrows and its lateral extent broadens. 
Finally, the opposite sides of the surface touch, with a 
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FIG. 2: Time-evolution of a slightly asymmetric neck cross- 
section. The initial distortion is a single n = 2 vibrational 
mode {A2/R0 = 0.004, Q2 = 7r/4). (a) The average radius 
rescaled by the initial average radius i^o, as a function of 
(t* —t)/t^, the amount of time remaining until the point im- 
plosion singularity at t* . (b) Neck cross-section at successive 
times. Shapes rescaled by the average radius R(t). (c) Shape 
evolution when n — 2 linear stability dynamics is extrapo- 
lated as a series expansion in 7^, V variables down to contact 
(Eqn. (8)). 



common tangent at the point of contact, thus severing 
itself into two side-by-side lobes. This is the typical out- 
come when an azimuthal asymmetry is present. It is also 
qualitatively consistent with the "puncture" type break- 
up (Fig. lb) observed in the experiment by Keim et al [6 
if one assumes that the satellite bubble at the center of 
the "puncture" is created by viscous drag created by the 
drainage of air from the thin gap in the final instants of 
contact. More complicated, but less frequent, types of 
self-intersection will be analyzed in a future study. 

Fig. 3 shows how Uc/{dR{t)/dt)^ the rescaled veloc- 
ity at the contact location, varies as the dimensionless 
gap width d/Ro goes to 0. Initially, while d is large, 
the interface accelerates inwards faster than the average 
radius decreases, or Uc/{dR{t)/dt) > 1. As d decreases, 
this rescaled velocity reaches a maximum and then slows. 
Eventually the rescaled contact speed settles to a con- 
stant value as d goes to 0. In Fig. 3, the simulation 
curves starting from 3 different initial conditions clearly 
terminate with different contact speeds, showing that Uc 
is chosen by the initial state. Note also that the sim- 
ulation typically yields Uc/ {dR{t) /dt) ^ 0(1), i.e. the 
cross-section shrinks about as fast as the two sides of the 
interface collides at the moment of contact. 

It is possible to give a simple interpretation of the qual- 
itative features found in these simulations. The contact 
dynamics found in a typical run can be thought of as 
an "extrapolation" of the amplitude-preserving vibration 
dynamics excited in the linear stability regime. Specifi- 
cally, we can rewrite the linear stability results for n = 2 
vibration (Eqn. (6,7)) as a Taylor series expansion in 7^, 
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FIG. 3: Velocity at the contact location, Uc/{dR/dt), as 
a function of the gap width d/Ro. Simulation results for 
A2/R0 = 0.004 (circles), 0.011 (crosses) and 0.067 (triangles) 
with the same Q2 = 7r/4 display the same qualitative behav- 
ior but terminate with different Uc values. The dot-dash line 
is prediction from Eqn. (9) with A2/R0 = 0.004 (circles). 



V (Eqn. (8,9)). 

n{w,t) = 



1 



wR{t) 

-1 
wR{t) 



1 



1 



A2 cos[MR)] \ ^ 

R{t) J 
A2sm[(t)2{R)]\ 1 



(8) 
(9) 



As A2 becomes comparable with R{t)^ Eqn. (8,9) repre- 
sents a particular scheme for approximating how the vi- 
brations transform when nonlinearity is significant. For 
our problem, this scheme successfully reproduces all the 
qualitative features found from the full simulation. This 
is evident in Fig. 2c, where the hypothetical shape evo- 
lution (Eqn. (8)) predicts that the cross-section osculate 
and reproduces the successive shape changes. Similarly, 
the time-evolution of the speed at the contact location 
predicted by Eqn. (9) show the same qualitative trend 
as the full simulation (circles), though the prediction 
does over-estimate the value of the final contact speed. 
Mathematically, this agreement says that 7^, V are "nat- 
ural" variables for tracking the break-up dynamics, since 
the time-evolution assumes such a simple form in 7^, 
V. Physically, this supports a simple scenario for how 
the initial distortion control the final contact dynamics: 
when the initial distortion is dominated by a single vibra- 
tional mode, the final contact is dominated by the same 
mode. Although the nonlinear interactions can and do 
generate new modes and cross-link different modes, in a 
typical situation neither mechanism gets enough time to 
cause qualitative changes. 

Fig. 4 plots Rc, the average radius of the cross- 
section at the moment of contact, as a function of the 
initial shape asymmetry A2/Ro- The prediction from 
Eqns. (8,9) (lines) reproduce both the contact orienta- 
tion and, roughly, the values of Rc found from the sim- 
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FIG. 4: Average radius Rc/Ro at the moment of contact as a 
function of the initial distortion amplitude A2/R0. The gap 
in the numerical result corresponds to initial conditions that 
evolve into more complicated forms of self intersection. 



ulations. Overall, Rc decreases with ^2, since smaller 
initial asymmetry allows the axisymmetric implosion to 
proceed further. The step-like structure associated with 
successive changes in the contact orientation reflects the 
fact that the final dynamics is dominated by the initial 
n = 2 distortion. Since a contact due to n = 2 vibration 
alone can only occur along two distinct orientations. 

It is worth noting that the finite-time contact is nei- 
ther a physical singularity, a divergence in the velocity or 
the pressure, nor a mathematical singularity, correspond- 
ing to one or more singularities of the mapping function 
crossing the unit circle on the complex plane [26l [28l [29] . 
Instead contact is a point where the mapping function be- 
comes multivalued. Mathematically, the bubble break-up 
problem occurs within a phase space which has a natural 
boundary, a "wall" separating singly-connected shapes 
from the multiply-connected ones. A contact corresponds 
to a time evolution that intersects a point on the "wall" 
in a finite amount of time. Because the time-evolution 
involves vibrations in the cross-section shape, both the 
speed and the average radius at contact has an intricate 
variation with respect to the initial amplitude and phase. 

In conclusion, we have conducted a numerical study of 
asymmetric bubble disconnection. The simulations show 
that an initial distortion excites vibrations which persist 
with essentially the same amplitude. The generic out- 
come is a finite-time contact. The distorted interface 
osculates with itself, forming a contact at a finite speed 
whose value depends on the initial condition. 
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